Viromes of Plants Determined by High-Throughput Sequencing of Virus-Derived siRNAs.

Plants growing in open airfields can be infected by several viruses even as a multiple infection. Virus infection in crops can lead to a serious damage to the harvest. In addition, virus presence in grapevine, fruit trees, and tuberous vegetables, propagated vegetatively affects the phytosanitary status of the propagation material (both the rootstock and the variety) having profound effect on the lifetime and health of the new plantations. The fast evolution of sequencing techniques provides a new opportunity for metagenomics-based viral diagnostics. Small interfering (si) RNAs produced by the RNA silencing-based host immune system during viral infection can be sequenced by high-throughput techniques and analyzed for the presence of viruses, revealing the presence of all known viral pathogens in the sample and therefore opening new avenues in virus diagnostics. This method is based on Illumina sequencing and bioinformatics analysis of virus-derived siRNAs in the host. Here we describe a protocol for this challenging technique step by step with notes, to ensure success for every user.

Clematis vitalba Is a Natural Host of the Novel Ilarvirus, Prunus Virus I.

Clematis vitalba L. is a climbing shrub and a pioneer plant in abandoned orchards or vineyards that are widespread in temperate climate zones. In past years, several viruses infecting the Clematis species have been identified, including different ilarviruses. Prunus virus I (PrVI) is a recently described ilarvirus, which has been shown to infect sweet cherries and peaches in Greece. Moreover, its presence has been detected in ornamental Clematis in Russia. In the present work, we analyzed the virome of wildly growing C. vitalba plants from Hungary, Slovakia and Croatia showing different kinds of symptoms using high-throughput sequencing (HTS) of small RNAs or ribodepleted RNAs. Applying HTS enabled us to identify the presence of PrVI in C. vitalba, and the bioinformatic analyses were further validated with RT-PCR using PrVI-specific primers and Sanger dideoxy sequencing. Nearly full genome sequences of all three viral RNAs of one Hungarian, two Slovak and one Croatian isolate were determined. Their phylogenetic analysis showed high similarity to each other and to other PrVI isolates described from Central Europe. As the sampled plants were co-infected with other viruses, it is not possible to determine a direct correlation between the infection with PrVI and the observed symptoms. Analyses of different Prunus species in stock collection showed infection of several peach and sweet cherry varieties in Hungary. Our results expand the knowledge on the natural host range of PrVI and highlight the necessity to evaluate alternative plant hosts (even non-Prunus) of PrVI and the role of the virus in the etiology of the potential diseases.

The First Virome of a Russian Vineyard.

Among other pathogens, more than 80 viruses infect grapevine. The aim of this work was to study the virome diversity of grapevine viruses and mycoviruses of a vineyard using high-throughput sequencing technologies. The grapevine virome was studied in symptomatic vines of the Rkatsiteli cultivar (V. vinifera) collected at the vineyards of the Krasnodar Krai in Russia. Ribosomal-depleted total RNA and isolated small RNAs were used for library preparation and high-throughput sequencing. Six grapevine-infecting viruses and two viroids were validated by RT-PCR and analyzed phylogenetically. We identified the presence of grapevine leafroll-associated virus 3, grapevine Pinot gris virus, grapevine virus T, grapevine rupestris stem-pitting-associated virus, grapevine fleck virus, and grapevine rupestris vein feathering virus, as well as two viroids, grapevine yellow speckle viroid 1 and hop stunt viroid. We also studied the mycovirome of the vineyard and identified nine viruses with single-stranded positive-sense RNA genomes: alternaria arborescens mitovirus 1, botrytis cinerea mitovirus 1, botrytis cinerea mitovirus 2, botrytis cinerea mitovirus 3, botrytis cinerea mitovirus 4, sclerotinia sclerotiorum mitovirus 3, botrytis cinerea hypovirus 1, grapevine-associated narnavirus 1, and botrytis virus F. In addition, we identified botrytis cinerea hypovirus 1 satellite-like RNA and two single-stranded negative-sense RNA viruses. This is the first study of grapevine mycoviruses in Russia. The obtained result will contribute to the development of biocontrol strategies in the future.

Analysis of Virus-Derived siRNAs in Strawberry Plants Co-Infected with Multiple Viruses and Their Genotypes.

Plants can be infected with multiple viruses. High-throughput sequencing tools have enabled numerous discoveries of multi-strain infections, when more than one viral strain or divergent genomic variant infects a single plant. Here, we investigated small interfering RNAs (siRNAs) in a single strawberry plant co-infected with several strains of strawberry mottle virus (SMoV), strawberry crinkle virus (SCV) and strawberry virus 1 (StrV-1). A range of plants infected with subsets of the initial viral species and strains that were obtained by aphid-mediated transmission were also evaluated. Using high-throughput sequencing, we characterized the small RNA fractions associated with different genotypes of these three viruses and determined small RNA hotspot regions in viral genomes. A comparison of virus-specific siRNA (vsiRNA) abundance with relative viral concentrations did not reveal any consistent agreement. Strawberry mottle virus strains exhibiting considerable variations in concentrations were found to be associated with comparable quantities of vsiRNAs. Additionally, by estimating the specificity of siRNAs to different viral strains, we observed that a substantial pool of vsiRNAs could target all SMoV strains, while strain-specific vsiRNAs predominantly targeted rhabdoviruses, SCV and StrV-1. This highlights the intricate nature and potential interference of the antiviral response within a single infected plant when multiple viruses are present.

Small RNA Profiling of Aster Yellows Phytoplasma-Infected Catharanthus roseus Plants Showing Different Symptoms.

Micropropagated Catharantus roseus plants infected with 'Candidatus Phytoplasma asteris' showed virescence symptoms, witches' broom symptoms, or became asymptomatic after their planting in pots. Nine plants were grouped into three categories according to these symptoms, which were then employed for investigation. The phytoplasma concentration, as determined by qPCR, correlated well with the severity of symptoms. To reveal the changes in the small RNA profiles in these plants, small RNA high-throughput sequencing (HTS) was carried out. The bioinformatics comparison of the micro (mi) RNA and small interfering (si) RNA profiles of the symptomatic and asymptomatic plants showed changes, which could be correlated to some of the observed symptoms. These results complement previous studies on phytoplasmas and serve as a starting point for small RNA-omic studies in phytoplasma research.

Detection of Apple Hammerhead Viroid, Apple Luteovirus 1 and Citrus Concave Gum-Associated Virus in Apple Propagation Materials and Orchards in the Czech Republic and Hungary.

Grafting cultivars onto rootstocks is a widely used practice by the apple industry predominantly aimed at faster fruit bearing. Using high-throughput sequencing, we revealed the presence of recently described viral agents, namely apple hammerhead viroid (AHVd), apple luteovirus 1 (ALV-1), and citrus concave gum-associated virus (CCGaV), in germplasm collections and production orchards in the Czech Republic and Hungary. The HTS results were validated with RT-(q)PCR, and Northern blotting techniques. To obtain further insight about the presence of these agents, RT-PCR based surveys were carried out and showed their widespread presence alone or in mixed infections. The pathogens were present both in production areas and in feral samples. In addition, rootstock-to-scion transmission of ALV-1 and CCGaV was confirmed using commercial rootstock materials. Phylogenetic relationships based on partial sequences of distinct variants were also investigated. Furthermore, the rosy apple aphid was found to be ALV-1-positive, suggesting that it might be a potential vector of the virus.

Grapevine Pinot Gris Virus Is Present in Different Non-Vitis Hosts.

Grapevine Pinot gris virus (GPGV) was described in Italy using a metagenomic approach: next-generation sequencing of the virus-derived small RNAs. Since that time, it has been reported all over the world. The presence of GPGV is associated with grapevine disease, but most of the time, the disease is asymptomatic. Although the host range of this virus has not been investigated, it has been found in the non-Vitis hosts, Silene latifolia and Chenopodium album. We investigated the presence of GPGV in grapevine and other plant species growing as weeds in the vineyard. Using RT-PCR, we identified GPGV in seven non-Vitis hosts: Ailanthus, Asclepias, Crataegus, Fraxinus, Rosa, Rubus, and Sambucus. In the case of Rosa and Rubus, this finding was supported by Northern blot detection of the virus. GPGV strains in non-Vitis hosts belong to the asymptomatic clade, and are clustered according to their original geographic locations. The presence of GPGV in species other than grapevine shows that besides well-known vector and propagating material-based infections, other possible entry sites for the virus can exist, which have to be taken into consideration when developing reliable regulation strategies.

Viromes of Hungarian Peach Trees Identified by High-Throughput Sequencing of Small RNAs.

Peach trees can be infected with viruses and viroids. As we do not have efficient plant protection methods against these pathogens, the prevention of infection is crucial. Fruit trees are maintained by vegetative propagation. Planting material such as certified mother trees and rootstocks should be free from viruses and viroids, and this status has to be regularly checked to prevent infections. We surveyed certified peach trees for the presence of viruses and viroids using small RNA high-throughput sequencing (HTS), an unbiased virus diagnostic method. The results of the bioinformatic analysis of HTS were validated by other molecular methods including RT-PCR, Northern blot hybridization and loop-mediated isothermal amplification (LAMP). We found the presence of plum pox virus and peach latent mosaic viroid (PLMVd) in the vector-free isolator houses, whose presence should be regularly tested. Moreover, we detected frequent infection with recently described viruses such as nectarine stem pitting-associated virus and peach-associated luteovirus (PaLV). During the survey, PLMVd and PaLV were detected for the first time in Hungary. The analysis of the presenting virus variants and possible sources of infection suggests that the source of the viral infection could be the infected propagating material. Our study emphasizes the importance of using sensitive and trustworthy diagnostic techniques to be able to detect viral infections and successfully prevent their spread by propagation material.

Controlled RISC loading efficiency of miR168 defined by miRNA duplex structure adjusts ARGONAUTE1 homeostasis.

Micro RNAs (miRNAs) are processed from precursor RNA molecules with precisely defined secondary stem-loop structures. ARGONAUTE1 (AGO1) is the main executor component of miRNA pathway and its expression is controlled via the auto-regulatory feedback loop activity of miR168 in plants. Previously we have shown that AGO1 loading of miR168 is strongly restricted leading to abundant cytoplasmic accumulation of AGO-unbound miR168. Here, we report, that intrinsic RNA secondary structure of MIR168a precursor not only defines the processing of miR168, but also precisely adjusts AGO1 loading efficiency determining the biologically active subset of miR168 pool. Our results show, that modification of miRNA duplex structure of MIR168a precursor fragment or expression from artificial precursors can alter the finely adjusted loading efficiency of miR168. In dcl1-9 mutant where, except for miR168, production of most miRNAs is severely reduced this mechanism ensures the elimination of unloaded AGO1 proteins via enhanced AGO1 loading of miR168. Based on this data, we propose a new competitive loading mechanism model for miR168 action: the miR168 surplus functions as a molecular buffer for controlled AGO1 loading continuously adjusting the amount of AGO1 protein in accordance with the changing size of the cellular miRNA pool.

High-Throughput Sequencing of Small RNAs for Diagnostics of Grapevine Viruses and Viroids in Russia.

The use of high-throughput sequencing (HTS) technology has led to significant progress in the identification of many viruses and their genetic variants. In this study, we used the HTS platform to sequence small RNAs (sRNAs) of grapevine to study the virome. Isolation of RNA was performed using symptomatic grapevines collected from commercial vineyards in Krasnodar Krai in 2017-2018. To determine the viromes of vineyards, we used an integrated approach that included a bioinformatic analysis of the results of sRNA HTS and the molecular method RT-PCR, which made it possible to identify 13 viruses and 4 viroids. Grapevine leafroll-associated virus 4 (GLRaV-4), Grapevine Syrah Virus-1 (GSyV-1), Raspberry bushy dwarf virus (RBDV), Australian grapevine viroid (AGVd), and Grapevine yellow speckle viroid 2 (GYSVd-2) were identified for the first time in Russia. Out of 38 samples analyzed, 37 had mixed infections with 4-11 viruses, indicating a high viral load. Analysis of the obtained sequences of fragments of virus genomes made it possible to identify recombination events in GLRaV-1, GLRaV-2, GLRaV-3, GLRaV-4, GVT, GPGV, GRSPaV, GVA, and GFLV. The obtained results indicate a wide spread of the viruses and a high genetic diversity in the vineyards of Krasnodar Krai and emphasize the urgent need to develop and implement long-term strategies for the control of viral grapevine diseases.

Variable Populations of Grapevine Virus T Are Present in Vineyards of Hungary.

Grapevine virus T (GVT) is a recently described foveavirus, which was identified from a transcriptome of a Teroldego grapevine cultivar in 2017. Recently, we surveyed vineyards and rootstock plantations in Hungary using small RNA (sRNA) high-throughput sequencing (HTS), at a time when GVT had not yet been described. A re-analysis of our sRNA HTS datasets and a survey of grapevines by RT-PCR revealed the presence of GVT in most of the vineyards tested, while at rootstock fields its presence was very rare. The presence and high variability of the virus in the country was confirmed by sequence analysis of strains originating from different vineyards. In this study, we demonstrate the presence of GVT in Hungary and show its high diversity, suggesting that GVT presence may not seriously affect grapevine health and that it could have been present in European vineyards for a long time as a latent infection.

HTS-Based Monitoring of the Efficiency of Somatic Embryogenesis and Meristem Cultures Used for Virus Elimination in Grapevine.

Meristem culture and somatic embryogenesis are effective tools for virus elimination of vegetatively propagated crops including grapevine (Vitis vinifera L.). While both have been shown to be useful to eliminate the main grapevine viruses, their efficiency differs depending on the virus and grapevine variety. In our work, we investigated the efficiency of these two virus elimination methods using small RNA high-throughput sequencing (HTS) and RT-PCR as virus diagnostics. Field grown mother plants of four clones representing three cultivars, infected with different viruses and viroids, were selected for elimination via somatic embryogenesis (SE) and meristem culture (ME). Our results show for the first time that using SE, elimination in mother plants was effective for all viruses, i.e., grapevine rupestris vein feathering virus (GRVFV), grapevine Syrah virus 1 (GSyV-1), Grapevine virus T (GVT) and grapevine Pinot gris virus (GPGV). This study also confirms previous studies showing that SE is a possible strategy for the elimination of GFkV, GRSPaV, HSVd, and GYSVd-1. Our results demonstrate that the efficacy of virus elimination via SE is relatively high while the purging of viroids is lower. Our work provides evidence that the efficiency of SE is comparable to that of the technically difficult ME technique, and that SE will offer a more effective strategy for the production of virus-free grapevine in the future.

Plasma activated water triggers plant defence responses.

Nowadays, one of the main challenges is moving towards an eco-sustainable agriculture, able to preserve the food production through a reduced use of pesticides. Current global food sustenance by intensive agriculture is mainly based on economic crop monocultures and drastically reduces the biodiversity, increasing the yield losses due to the presence of biotic and abiotic stresses. A technology based on plasma activated water (PAW), characterized by the presence in liquid of reactive oxygen and nitrogen species, was tested to try to ensure yield stability also enhancing the plant resistance responses and to promote an eco-sustainable management of plant diseases. In PAW-treated micropropagated periwinkle shoots, periwinkle and grapevine plants, qRT-PCR and small RNAs high-throughput sequencing were used to analyse the differential expression of genes involved in the major plant defence pathways. The results indicate that PAW treatment enhances the plant defence responses and provide an encouraging framework for future applications in plant disease management programs.

Millet Could Be both a Weed and Serve as a Virus Reservoir in Crop Fields.

Millet is a dangerous weed in crop fields. A lack of seed dormancy helps it to spread easily and be present in maize, wheat, and other crop fields. Our previous report revealed the possibility that millet can also play a role as a virus reservoir. In that study, we focused on visual symptoms and detected the presence of several viruses in millet using serological methods, which can only detect the presence of the investigated pathogen. In this current work, we used small RNA high-throughput sequencing as an unbiased virus diagnostic method to uncover presenting viruses in randomly sampled millet grown as a volunteer weed in two maize fields, showing stunting, chlorosis, and striped leaves. Our results confirmed the widespread presence of wheat streak mosaic virus at both locations. Moreover, barley yellow striate mosaic virus and barley virus G, neither of which had been previously described in Hungary, were also identified. As these viruses can cause severe diseases in wheat and other cereals, their presence in a weed implies a potential infection risk. Our study indicates that the presence of millet in fields requires special control to prevent the emergence of new viral diseases in crop fields.

Local Aphid Species Infestation on Invasive Weeds Affects Virus Infection of Nearest Crops Under Different Management Systems - A Preliminary Study.

In the present study, we conducted field surveys to detect the population density of the most important invasive weed species and their associated virus vectoring aphids in crops grown under high input field (HIF) vs. low-input field (LIF) conditions, with and without fertilizers and pesticides. The most frequent invasive weed species were annual fleabane, Erigeron annua (L.), Canadian horseweed, Erigeron canadensis (L.) and Canadian goldenrod, Solidago canadensis (L.). These species were predominantly hosts of the aphids Brachycaudus helichrysi and Aulacorthum solani under both management systems. The 13% higher coverage of E. annua under LIF conditions resulted in a 30% higher B. helichrysi abundance and ∼85% higher A. solani abundance compared with HIF conditions. To reveal the incidence of virus infection in crop plants and invasive weeds, high-throughput sequencing of small RNAs was performed. Bioinformatics analysis combined with independent validation methods revealed the presence of six viruses, but with strikingly different patterns under LIF and HIF conditions. Their presence without symptoms in invasive weeds and crop plants supports the necessity of employing new approaches to those currently employed in invasive weed management. These findings also suggest that invasive weeds could serve as hosts for local aphid species and reservoirs for plant pathogenic viruses, both under low and high input management systems. In this light, as here demonstrated, viruses transmitted by local aphid species were found to differ between the management systems; hence, the importance of B. helichrysi and A. solani as virus vectors in particular clearly needs to be re-evaluated. Altogether, we accept that the present study is a pilot one and individual virus vectoring of aphids still needs to be directly tested. Even so, it represents one of the first contributions to this particular area, and thereby paves the way for further similar applied research in the future.

Complete Sequence, Genome Organization and Molecular Detection of Grapevine Line Pattern Virus, a New Putative Anulavirus Infecting Grapevine.

: Grapevine line pattern virus (GLPV) was first described 30 years ago in Hungary. The lack of its genomic sequences and of an available antiserum made its detection impossible in other parts of the world. Three different high-throughput sequencing (HTS) protocols applied on a GLPV-infected vine allowed the construction of the full genome sequence of this virus. It includes three RNA segments, encoding four proteins: methyltransferase-helicase (1a), RNA-dependent RNA polymerase (2a), movement protein (3a) and coat protein (3b). The obtained sequences were used to design specific primers for its detection by RT-PCR and Northern blot hybridization, respectively. These diagnostic methods were used to test the presence of GLPV in graft-inoculated plants and in 220 grapevine accessions of different Mediterranean origins. The three RNAs-encoding proteins of GLPV shared a very high amino acid identity with those of hop yellow virus, a tentative member of the Anulavirus genus, leaving no doubt that both are two isolates of the same viral species. A circular RNA originating from the RNA2 was found, for which an alternative silencing suppressor role is hypothesized. Further investigation is needed to determine this possibility and also the host range and pathological significance of the virus.

Analysis of a novel RNA virus in a wild northern white-breasted hedgehog (Erinaceus roumanicus).

Tombusviruses are generally considered plant viruses. A novel tombus-/carmotetravirus-like RNA virus was identified in a faecal sample and blood and muscle tissues from a wild northern white-breasted hedgehog (Erinaceus roumanicus). The complete genome of the virus, called H14-hedgehog/2015/HUN (GenBank accession number MN044446), is 4,118 nucleotides in length with a readthrough stop codon of type/group 1 in ORF1 and lacks a poly(A) tract at the 3' end. The predicted ORF1-RT (RdRp) and the capsid proteins had low (31-33%) amino acid sequence identity to unclassified tombus-/noda-like viruses (Hubei tombus-like virus 12 and Beihai noda-like virus 10), respectively, discovered recently in invertebrate animals. An in vivo experimental plant inoculation study showed that an in vitro-transcribed H14-hedgehog/2015/HUN viral RNA did not replicate in Nicotiana benthamiana, Chenopodium quinoa, or Chenopodium murale, the most susceptible hosts for plant-origin tombusviruses.

Differential gene expression and physiological changes during acute or persistent plant virus interactions may contribute to viral symptom differences.

Viruses have different strategies for infecting their hosts. Fast and acute infections result in the development of severe symptoms and may cause the death of the plant. By contrast, in a persistent interaction, the virus can survive within its host for a long time, inducing only mild symptoms. In this study, we investigated the gene expression changes induced in CymRSV-, crTMV-, and TCV-infected Nicotiana benthamiana and in PVX- and TMV-U1-infected Solanum lycopersicum plants after the systemic spread of the virus by two different high-throughput methods: microarray hybridization or RNA sequencing. Using these techniques, we were able to clearly differentiate between acute and persistent infections. We validated the gene expression changes of selected genes by Northern blot hybridization or by qRT-PCR. We show that, in contrast to persistent infections, the drastic shut-off of housekeeping genes, downregulation of photosynthesis-related transcripts and induction of stress genes are specific outcomes with acute infections. We also show that these changes are not a consequence of host necrosis or the presence of a viral silencing suppressor. Thermal imaging data and chlorophyll fluorescence measurements correlated very well with the molecular changes. We believe that the molecular and physiological changes detected during acute infections mostly contribute to virus symptom development. The observed characteristic physiological changes associated with economically more dangerous acute infections could serve as a basis for the elaboration of remote monitoring systems suitable for detecting developing virus infections in crops. Moreover, as molecular and physiological changes are characteristics of different types of virus lifestyles, this knowledge can support risk assessments of recently described novel viruses.

Virus Detection by High-Throughput Sequencing of Small RNAs: Large-Scale Performance Testing of Sequence Analysis Strategies.

Recent developments in high-throughput sequencing (HTS), also called next-generation sequencing (NGS), technologies and bioinformatics have drastically changed research on viral pathogens and spurred growing interest in the field of virus diagnostics. However, the reliability of HTS-based virus detection protocols must be evaluated before adopting them for diagnostics. Many different bioinformatics algorithms aimed at detecting viruses in HTS data have been reported but little attention has been paid thus far to their sensitivity and reliability for diagnostic purposes. Therefore, we compared the ability of 21 plant virology laboratories, each employing a different bioinformatics pipeline, to detect 12 plant viruses through a double-blind large-scale performance test using 10 datasets of 21- to 24-nucleotide small RNA (sRNA) sequences from three different infected plants. The sensitivity of virus detection ranged between 35 and 100% among participants, with a marked negative effect when sequence depth decreased. The false-positive detection rate was very low and mainly related to the identification of host genome-integrated viral sequences or misinterpretation of the results. Reproducibility was high (91.6%). This work revealed the key influence of bioinformatics strategies for the sensitive detection of viruses in HTS sRNA datasets and, more specifically (i) the difficulty in detecting viral agents when they are novel or their sRNA abundance is low, (ii) the influence of key parameters at both assembly and annotation steps, (iii) the importance of completeness of reference sequence databases, and (iv) the significant level of scientific expertise needed when interpreting pipeline results. Overall, this work underlines key parameters and proposes recommendations for reliable sRNA-based detection of known and unknown viruses.

Small RNA NGS Revealed the Presence of Cherry Virus A and Little Cherry Virus 1 on Apricots in Hungary.

Fruit trees, such as apricot trees, are constantly exposed to the attack of viruses. As they are propagated in a vegetative way, this risk is present not only in the field, where they remain for decades, but also during their propagation. Metagenomic diagnostic methods, based on next generation sequencing (NGS), offer unique possibilities to reveal all the present pathogens in the investigated sample. Using NGS of small RNAs, a special field of these techniques, we tested leaf samples of different varieties of apricot originating from an isolator house or open field stock nursery. As a result, we identified Cherry virus A (CVA) and little cherry virus 1 (LChV-1) for the first time in Hungary. The NGS results were validated by RT-PCR and also by Northern blot in the case of CVA. Cloned and Sanger sequenced viral-specific PCR products enabled us to investigate their phylogenetic relationships. However, since these pathogens have not been described in our country before, their role in symptom development and modification during co-infection with other viruses requires further investigation.